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Abstract. We have searched for Extremely Red Objects (EROs) around faint mid-IR selected galaxies in ELAIS 
fields. We find a significant overdensity, by factors of 2 to 5, of these EROs compared to field EROs in the same 
region and literature random field ERO counts. The excess is similar to that found previously in the fields of 
known high redshift quasars and AGN. While with the present data it cannot be definitely shown whether the 
overdensity is physically connected to the mid-IR source, a derived radial distribution does suggest this. The 
fraction of EROs among Tf-selected galaxies is high in the mid-IR fields in agreement with the picture that the 
EROs responsible for the overdensity are members of high redshift clusters of galaxies. We find R — K > 5 
selected EROs to be more clustered around the mid-IR galaxies than I — K > A EROs, though statistics are 
weak because of small numbers. However, this would be consistent with a cluster/galaxy group scenario if, as 
we argue, the former colour selection finds preferentially more strongly clustered early type galaxies, whereas the 
latter selection includes a larger fraction of dusty EROs. Finally, using the mid-IR data, we are able to limit the 
fraction of ULIRG type very dusty EROs at A' < 18 magnitude to less than ~ 10% of the total ERO population. 

Key words, galaxies: evolution - infrared: galaxies - cosmology: observations - galaxies: clusters: general 



1. Introduction 

Galaxy formation and evolution is undoubtedly one of the 
major issues under debate in modern astronomy and cos- 
mology. On the observational side, there has been a great 
effort to gather samples of galaxies against which the suc- 
cesses and short-comings of galaxy formation models can 
be weighed. Two useful approaches are determining the 
space densities and ages of the most massive galaxies 
at and beyond z 1, and identifying clusters or proto- 
clusters of galaxies at those redshifts. 

The definition of Extremely Red Objects (EROs) for 
example by colours of i? — iiT > 5, / — i^T > 4, is motivated 
by this selection pin-pointing the crucial evolved and most 
massive galaxies at z > 1 from optical and near-infrared 
extragalactic imaging surveys. The models are often pre- 
sented as two extremes, where on one side there are mono- 
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lithic collapse models (e.g. Eggen, Lynden-Bell & Sandage 
119621 Larson [T975|l with pure luminosity evolution, and 
on the other hierarchical assembly scenarios (e.g. White 
& ReesETHI White & Frenk lTMTl Somerville & Primack 
Cole et al.|inD(T). The predicted numbers of evolved 
galaxies at high redshift have been very different, at least 
until recently (see eg. Somerville et al. 2004), and thus 
in principle easy to test against observations. However, in 
addition to the fact that reality possibly is some mix of 
models (e.g. Roche et al. l2002|l . the results from observed 
numbers of high-z ellipticals have also been inconclusive. 

There are two difficulties in interpreting the ERO num- 
ber count results. First of all, the extreme red colours 
can be attributed also to dust reddened starburst galax- 
ies thus contaminating the old distant early type sam- 
ples (see e.g. recent papers by Cimatti et al. 120031 Yan 
& Thompson .2003.,: Takata et al. 120031 and references 
therein). According to Cimatti et al. H2002|l and Small 
et al. 30 - 60% of EROs at typical ERO bright- 

nesses of if ~ 18 — 20 might be dusty. At brighter levels, 
K < 17.5, we find a smaller percentage < 20% of dusty 
EROs (Vaisanen & Johansson 2004, hereafter Paper I). 
Recent HST morphological studies have complicated the 
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picture somewhat by finding large fractions of disk galaxy 
EROs (Yan & Thompson HUU^ Moustakas et al. [SUMl 
Gilbank et al. l2(JU3fl . In addition to morphological, spec- 
troscopic, or colour property methods, the separation of 
dusty and early type EROs can also be done using de- 
tections at longer wavelengths, where ellipticals are not 
expected to radiate much. For attempts in sub-mm to 
radio see e.g. Mohan et al. (PPH^ . Small et al. (P02), 
Wehner, Barger, & Kneib pIT)^ . Webb et al. ("5004). In 
Paper I we used for the first time a mid-IR data-set to 
select dusty EROs - the difference of ellipticals and star- 
forming galaxies is very large also in the mid-IR, and e.g. 
the identification problems of large sub-mm beams can be 
avoided. However, the exact relative fractions of different 
types of EROs as a function of magnitude is still an open 
question. 

The second reason for varying numbers of EROs is 
their strong clustering (Daddi et al. I2000al Firth et al. 
nnn^ Roche et al. 12002)) . which is in fact expected if EROs 
are the counterparts of present day massive ellipticals. 
This clustering ties the ERO studies together with the 
second approach to understanding galaxy and large scale 
structure formation, namely the search for high redshift 
clusters of galaxies (see e.g. Olsen et al.'lSQO: Gonzalez et 
al. Onrm BahcaU et aL lMini Mullis et al. ,2003,) . It is fea- 
sible that EROs are found in clusters and concentrations 
of EROs can be used to locate high-z clusters of galaxies. 

Significant overdensities of EROs, compared to large 
field surveys, have been found in high-z AGN and QSO 
fields (e.g. McCarthy, Persson, & West HMS, Hu & 
R.idgewav ll994l Dey, Spinrad, & Dickinson ll995t Cimatti 
et al.Ennni chapman, McCarthy & Persson 1211101 Hall et 
al. l^nUn Sanchez & Gonzalez-Serrano TMM Best et al. 
l^nnSl Wold et al. HUU^ . As discussed in these studies, 
it is not yet clear whether the overdensity is necessar- 
ily physically connected to the AGN: it is also possible 
that the EROs would be part of foreground overdensities 
of galaxies which enhance the probability of detection of 
these AGN via gravitational lensing (see e.g. discussion in 
Wold et al. ,2003,) . Nevertheless, in either case the EROs 
can be used to identify candidates for high redshift galaxy 
clusters, and thus once again present themselves as an im- 
portant window into the galaxy and structure formation 
era of the Universe. 

We present here a NIR survey performed in the 
European Large Area ISO Survey (ELAIS; Oliver et al. 
I2(}()()|l fields, and derive surface densities of EROs. Rather 
than a survey with large continuous sky coverage (as in 
Paper I) the deeper NIR data-set of the present paper con- 
sists of numerous individual pointings around faint ELAIS 
ISO CAM detections or blank fields. We make use of 
the newly available Final band merged ELAIS Catalogue 
(Rowan-Robinson et al. 120041 Vaccari et al. 2004) and a 
publicly available optical survey. 

In particular, as the main aim of this paper, we com- 
pare the distribution of EROs in "random" fields to those 
fields which are occupied by faint mid-IR detections. We 
are thus able to determine whether there are similar over- 



density biases of EROs as in the case of verified high-z 
AGN. Similar results are expected at least to some ex- 
tent, since a significant fraction of faint ISOCAM sources 
appear to be AGN at z > 1 (e.g. Fadda et al. 20021 La 
Franca et al.[2004, Johansson, Vaisanen & Vaccari 2004, 
in prep.). As mentioned, locations of ERO overdensities 
may serve as candidates for high-redshift clusters and 
groups of galaxies, whose history in turn place strong con- 
straints on structure formation scenarios. A demonstrated 
link between ERO concentrations and specific classes of 
sources (in this case AGN or certain mid-IR galaxies) 
would then further facilitate finding the crucial clusters 
using data such as large MIR surveys with the Spitzer 
Space Telescope. 

The structure of this paper is as follows. In Sect. 2 we 
present the observational data. In Sect. 3 the properties 
of the ERO samples are defined and analysed, and Sect. 4 
discusses the overdensity of EROs. We assume throughout 
this paper a flat {Qq ~ 1) cosmology with = 0.3, 
nA = 0.7 and Ho = 70 kms-^Mpc-i. 

2. Observations 

2.1. ISO data 

The ISO data used in this study comes from the Final 
band-merged ELAIS catalogue (Rowan-Robinson et al. 
2004, Vaccari et al. 2004). For further description of the 
ELAIS project, observations, data reduction, source ex- 
traction and extragalactic source counts see Oliver et al. 
tiOOO), Serjeant et al. (2000), and Efstathiou et al. pWill . 
The observations used in this study were primarily made 
using the ISOCAM (Cesarsky et al. [TUU^ LW3 filter at 
15 fim. 

2.2. Near-Infrared data 

The observations presented here were carried out with 
the 3.0-m NASA Infrared Telescope Facihty (IRTF) on 
Mauna Kea, and were conducted in the ELAIS fields Nl 
and N2, centered at (a, (5) = 16h 09m 00s, 54° 40' 00" and 
{a,S) = 16h 36m 00s, 41° 06' 00" J2000.0 respectively The 
data were acquired during five photometric nights in June 
1998, and the seeing varied between 0.7" — 1.0". The 
data set consists of images taken at (very preliminary) 
ISOCAM detection positions where no obvious counter- 
parts were found in the Digital Sky Survey (DSS). The 
sample used here consists of a total of 33 images. Of these, 
9 were taken from the ELAIS Nl region and 24 from the 
ELAIS N2 region. The total surveyed area is 73.6 arcmin^. 
These areas are small non-continuous regions of the entire 
ELAIS Nl and N2 fields, which cover a continuous area of 
> 2 deg2. 

We used the NFSCam instrument containing a 256 x 
256 format InSb detector array. A 0.3" pixel scale was 
used for all observations. Images were constructed from co- 
added dithered frames using standard NIR data reduction 
methods (such as in DIMSUM) written as IDL procedures. 
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The final images have a field of view of 1.5' x 1.5'. The total 
integration times were in the range of 20 - 30 minutes and 
the limiting magnitude for the deepest images was around 
K = 20.0. 

2.3. Optical data 

Optical photometric data were obtained from the Isaac 
Newton Telescope (INT) Wide Field Survey (WFS). 
The WFS data are publicly available on the Cambridge 
Astronomical Survey Unit (CASU) homepage^. See 
McMahon et al. for a review of the INT Wide Field 

Survey Project. We use data consisting of Vega-based 
g' ,r' ,i' , Z band photometry (the apostrophes are dropped 
henceforth). The filters are similar to the SDSS (Sloan 
Digital Sky Survey) filters (Fukugita et al. '1996). We 
heavily rely on the r, i bands, and use the following filter 
transformations to the standard Cousins system (Landolt 
WFS webpage^): r-R = 0.275(i? - /) + 0.008 and 
i — I = 0.211(i? — /). Our adopted 5ct detection limits for 
a 1" seeing are g = 25.0, r = 24.1, i = 23.2 and z = 22.0 
(WFS webpage). For further details on the pipeline pro- 
cessing of INT wide field survey data consult Irwin & 
Lewis H2001fl and Gonzalez-Solares et al. H2004|l for dis- 
cussion of the WFS data in ELAIS fields in particular. 

2.4. Other data 

In addition to 90 and 175 fim, and 20 cm data available 
in the ELAIS catalogue, we checked several other separate 
surveys in the same region. A region of the ELAIS N2 field 
has been surveyed in the sub-millimetre with SCUBA to a 
3.5ct limit F{850nm) ~ 8 mJy by Scott et al. ij^UU^ . In the 
x-ray, the ROSAT (0.1-2 keV) all-sky survey was cross- 
correlated with ELAIS sources (Basilakos et al. l2002|l and 
deep 75 ks Chandra observations were performed in Nl 
and N2 regions (Almaini et al. :2003: Manners et al.[2D03). 
None of these surveys produced overlapping sources with 
the final list of EROs analyzed in this paper. In the same 
N2 area Roche et al. H2()()2|l have performed a deeper study 
of the clustering, number counts and morphology of EROs. 
Only one of our IRTF images overlaps with this very well 
studied region. 

2.5. Astrometry and photometry 

The final astrometry of the IRTF data was obtained by 
matching the visible sources with with Hubble Guide Star 
Catalog 2 objects (GSC2, available at the GSC home- 
page'^) giving us an accuracy of « 0.3". The INT astrom- 
etry was derived using GSC objects, which results in an 
external astrometric accuracy of . 5" — 1 . 0" (Irwin & Lewis 

mm . 



^ http://archive.ast.cam.ac.uk/ 

^ http:/ /www. ast.cam.ac.uk/~wfcsur/index.php 

^ http:/ /www-gsss. stsci.edu/gsc/GSChome. htm 



The photometry and calibration of the optical INT 
Wide Field Survey data were obtained from the archives 
(Irwin & Lewis l^ron Gonzalez-Solares et al OUMIl . 

All our near-IR photometry is performed using the 
SExtractor software (v. 2. 2.1 and v. 2. 3; Berlin & Arnouts 
11996(1 . For detecting objects in the IRTF data we used the 
following SExtractor settings: a 1.25f7sky (translating to 
K ^ 20.7 mag arcsec^^ for integration times > 20 min and 
K ^ 20.4 mag arcsec~^for those with 10-15 min) detec- 
tion threshold in a minimum area of 5 pixels (0.45 arcsec^), 
with a detection filter of 2.5 pixels Gaussian FWHM, ap- 
proximately matching the seeing. 

For total magnitudes, the SExtractor BEST- 
magnitude was used - it is most often the Kron-magnitude 
fKron [T?)8()|l . where an elliptical aperture for photometry 
is defined by the shape and size of the detected object. 
We found these magnitudes to be robust over a wide 
range of magnitudes and source profiles and contain 
(nearly) all of the object flux. In very crowded regions the 
BEST-magnitude can also be the isophotal magnitude 
instead. However, since the WFS data is given in 2.4" 
diameter apertures, we calculated corresponding aperture 
magnitudes from the IRTF images. All optical-near-IR 
colours given in this paper are thus calculated using 
matched apertures. 

As mentioned, the deepest parts of the NIR data reach 
K — 20. Noticeable incompleteness starts affecting NIR 
source counts at K > 19 mag. Since we are interested in 
the very reddest sources, and the detection limits of the 
optical bands are r w 24.1, i w 23.2 we will not consider 
NIR sources fainter than if = 18.7 in this paper. We thus 
do not need any completeness corrections to number count 
data. 

The magnitude calibration was determined from 
short exposures of UKIRT faint standard stars and 
HST/NICMOS IR standards. From these observations a 
magnitude zeropoint of 21.47 was derived, which is consis- 
tent with the IRTF/NFSCAM manual value. We estimate 
our K-hand photometry to be accurate within 0.06 mag. 
We also note that the ERO list in Roche et al. H2002|l has 
4 objects within one of our IRTF frames. Of these, 3 have 
fully consistent magnitudes though the fourth is 0.7 mag 
brighter in our catalogue. 

3. Construction of the ERO samples 

3.1. Definition of EROs 

Numerous different selection criteria have been defined for 
EROs, including R-K>6, R~K> 5.3, R - K > 5, 
I — K > 4 with if-magnitude upper limits from 18 to 20. 
In this paper we use the following definition for EROs: 
r — K > 5.5 and/or i — K > 4.4. Using the colour trans- 
formations given in Sec. 12.31 above, and typical colours of 
our objects, i? — J « 1.4 — 2 (see e.g. Fig. 2 of Paper I), 
it is seen that our choices correspond to the widely used 
R-K>5, I-K>4 ERO selection criteria. AU ERO 
selections are designed for selecting early type galaxies at 
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z > 1 - however, we note that since the model colours 
that different authors use vary, there might easily be 0.1- 
0.3 mag differences in the colours at z ~ 1. The elliptical 
galaxy models we used (GRASIL library; Silva et al. ll998|l 
become EROs at z w 1.1 — 1.2 with the criteria defined 
above. 

FigureHshows r — K and i — K model colours of several 
representative galaxies against redshift, along with the 
ERO criteria. Model SEDs are adopted from the GRASIL 
library (Silva et al. ri9i?5l a sample set of model SEDs 
and procedures to calculate more are available online^). 
Ordinary spirals never reach the red colours of EROs, 
while both ellipticals and reddened starbursts become 
EROs when seen beyond z ~ 1. For comparison, the colour 
of the prototype dusty ERO HRIO is also plotted as a 
function of redshift. The lowest panel shows the flux ratio 
fi5fim/ f2.2fim " the degeneracy in the red colours of old 
ellipticals and dusty starbursts is clearly broken. 

3.2. Matching of optical and NIR data 

The optical counterparts of ncar-IR sources were extracted 
from the INT WFS catalogue using a 1.5" search radius. 
Since the WFS catalogue contains multiple detections of 
the same source (due to overlapping CCD frames), optical 
counterparts within 0.5" were averaged. After the purg- 
ing, the brightest optical source of any remaining multi- 
ple counterparts was selected. There were 118/722 = 16% 
NIR sources which were undetected in the WFS catalogue. 
This is expected, since the IRTF survey goes somewhat 
deeper than the WFS data. For undetected sources the 
respective 5cr limiting magnitude was assigned as the lim- 
iting brightness. 

3.3. Star vs. galaxy separation 

Next, stars were separated from galaxies. Stellarity indices 
are available from both the NIR catalogue (the SExtractor 
CLASS parameter) and the WFS catalogue, where a flag 
defines galaxies, definite stars, and various degrees of un- 
certain stellarities. Since the IRTF data is deeper the NIR 
data provide more reliable classifications than the WFS. 

For the EROs we are interested in, the reddest and 
faintest objects in the lists, we found that a colour selec- 
tion works at least as well, if not better, than the morpho- 
logical separation. Please refer to Sect. 3.3 and Fig. 2 in 
Paper I, where a separating line of r—K > 2.16(r— i)-|-1.35 
for galaxies is derived. At the r — K > 5 regime there is 
no overlap between the populations. In essence, very late 
type cool stars (in particular of spectral type L) can have 
optical to NIR colours very similar to extragalactic EROs 
(e.g. Dahn et al. l2(JU2|l . However, the r — i colour is > 2 for 
these stars, whereas galaxies are bluer. Using this colour, 
extremely red stars can be discarded from the ERO sam- 
ples. 



* http:/ /web. pd.astro.it/granato/ 
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Fig. 1. Observed r — K and i — K and MIR/NIR colours 
as a function of redshift. Representative model galaxy 
SEDs are taken from GRASIL models (Silva et al. [T!?M| : 
an evolving Sb spiral (dotted line), an evolving ellipti- 
cal (solid curve), non-evolving prototypical starburst M82 
(dash-dot), and HRIO ERO (dashed). The dash-triple- 
dot curve in the lowest panel corresponds to the ULIRG 
ARP220. The horizontal lines show the colour definition 
of EROs in this paper. 



In practice, we used a combination of methods: At 
K < 18 we define the stellar sources according to NIR 
morphology (those with CLASS>0.5). At if > 18 those 
with CLASS>0.5 and a stellar colour are assigned as stars 
(if detected optically, only morphology used if no optical 
detection was available). We find 128 stars at K < 18 mag, 
which is in very good agreement with the prediction of 115 
stars for these fields from the Galactic foreground point 
source model SKY (Cohen et al. I19H4I see also Vaisanen 
et al. I2()0()|l . As for the colours of stars, approximately 
55% and 35% of objects redder than r — K > 5.5 (i.e. 
R — K > 5 are found to be stars at the limits if < 18 
and K < 18.5, respectively. These are very consistent 
with respective values calculated from the stellar distri- 
bution model of Jarrett fl992' Wold, private communica- 
tion; Wold et al. i2003jl . In summary, it should be noted 
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Fig. 2. An example of one of our IRTF i^T-band fields; the 
r-band data is from the WFS archive. There are many red 
sources in this field around ELAISC15_J163541. 7+405900 
(indicated by the arrow), and those satisfying our ERO 
criteria are circled. The four sources correspond to EROs 
#7, 6, 34, and 23 in TableEl starting from the lower left. 
The extremely red ERO-6 lies within 3a of positional un- 
certainty of the MIR source (see Section ing . 



that in contrast to fainter K > 19 ERO surveys, brighter 
ERO surveys are significantly contaminated by very cool 
stars in our own Galaxy. 

3.4. Extraction of EROs 

We searched for EROs from the resulting galaxy cata- 
logue according to the colour definitions given above in 
Section im Only NIR detections at 5a level and over were 
considered. As already mentioned, there were 118 sources 
which were not detected to the WFS survey limits: These 
are assigned magnitudes of r > 24.1 and i > 23.2, and thus 
may enter our ERO catalogue with colours r — K > 5.5 
and i — K > 4.4. We then went through the resulting list 
checking our NIR maps as well as the WFS CCD images 
for any obvious spurious objects. Ultimately there were 
27 EROs in the matched IRTF catalogue using the r — K 
criterion and 27 using the i — K limit — 14 EROs are 
common to lists resulting from both selection criteria, i.e. 
there were 40 EROs found in total (TableHJ. Note in par- 
ticular that 19/40 of the EROs have only upper limits in 
optical bands, ie. come from among those 118 sources de- 
tected only in the NIR. An example of an ERO field is 
shown in Fig. [3 

The total numbers of verified EROs are summarized 
in Tabled along with some other survey characteristics. 
Colour-magnitude plots are shown in Fig. |3| We also note 
that of the 4 common EROs with the Roche et al. H2002|l 
catalogue, one is excluded from our survey as a star and 
another because of detection levels (i.e. it is not detected in 
the WFS, but a. K > 18.7 magnitude does not necessarily 
make it an ERO as far as our data is concerned). 



4. Discussion 

4.1. Separation of dusty and early type EROs 

One of the motivations of this study is to attempt the 
se-paration of dusty EROs from the evolved early type 
population using mid-IR data. As seen in the lowest panel 
of Fig. n the /i5^im//2.2/im Aux ratio is below unity for 
an elliptical galaxy in the redshift range of z = — 3, 
whereas the ratio is ~ 10 for a typical dusty ERO. The 
fibfmi/ f2.2fmi ratio is thus in principle a very powerful 
discriminator between dustfree and dusty EROs. 

What fraction of the EROs, to a given flux limit, are 
detected in the mid-IR? Because our survey is not a pure 
field survey, but rather a targeted one, we cannot answer 
the question in general. Various limits can be obtained, 
however. 

One clear ISO-ERO, and two ambiguous cases were 
detected among the total of 40 IRTF EROs. In one of the 
ambiguous cases, an ISOCAM source lies within 3" of both 
an ERO and a spectroscopically confirmed non-ERO Syl 
galaxy at z — 1.150 which is selected as the counterpart 
in the ELAIS catalogue (Rowan-Robinson et al.P004). In 
the other ambiguous case a MIR source is associated with 
a. z — 0.1882 optical spectroscopic detection of a "spiral" 
galaxy in the ELAIS catalogue. However, the optical-NIR 
magnitudes, colours and size of this galaxy are that of 
Im type, and the strong 6.7 and 15 ^m emission detected 
would be difficult to understand (unless the source is a 
strong starburst or AGN, which should have been obvious 
from the spectrum). Instead, there is a very red ERO 6" 
away to which the IR SED would fit much better (the 
la positional uncertainty of the ELAIS MIR catalogue is 
w 2"). In fact, the preliminary ELAIS source in the source 
list used at the time of our NIR observations, was closer 
to this ERO than the z = 0.188 dwarf. 

If we assumed that the IRTF frames are a non-biased 
sample of the NIR sky, we would derive a lower limit of 
3 — 8% for the dusty EROs in the total ERO population 
(ignoring detection limits for the moment). In fact, the 
assumption is not totally incorrect, since only 1/3 of our 
frames contain a ISOCAM source (the reason being the 
prelimary source list). 

However, it is important to stress that the ELAIS mid- 
IR detection limit of 0.7 mJy at 15/im results in severe 
limitations. This detection limit corresponds to, for ex- 
ample, a starburst galaxy of K = 17.4 brightness, with a 
SED shape of M82 located at z = 1. This would already 
be very close to the ULIRG limit of Lju = 10^^ Lq (see 
also e.g. Elbaz et al. 12002(1 . and all our EROs are fainter 
than this in the NIR. As discussed in Paper I, deeper mid- 
IR imaging with e.g. the Spitzer/SIRTF mission, is clearly 
needed to detect counterparts of those dusty EROs which 
are not (ultra)luminous IR galaxies. 

As for the more extreme dusty sources, an ARP220 
SED-type object at z w 1 and K < 18 would have > 
1.3 mJy at 15/xm, and thus should be readily detected 
by the ELAIS survey. For comparison with other dusty 



6 



Vaisanen & Johansson: ERO overdensity around mid-IR galaxies 




Fig. 3. The K-hand magnitude vs. the r — K and i — K colours. Only galaxies are plotted. Verified EROs are overlaid 
with solid circles, and ISOCAM detections as solid triangles The objects lying on the r ~ 24.1 and i — 23.2 lines are 
sources which were not detected in the WFS and for which the Scr detection limit is assigned. Note that the total 
ivT-magnitude is used on the x-axis, whereas the colour is calculated with matching small apertures. 





N 


iV* 


Area arcmin^ 




<r- K > 


<i-K> 




EROs 


ERO/arcmin^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) (10) 


(11) (12) 


IRTF 


722 


128 


74 


19.5 


3.80 


2.90 


27 


27 14 


0.36 0.122-0.189 



Table 1. Columns (2) to (5) give the total numbers of objects and stellar sources in our survey, the covered area, and 
the characteristic limiting magnitude. The rest of the columns refer to galaxies only: columns (6) to (7) give mean 
colours, though note that these are lower limits since optical non-detections we ignored. Column (8) gives the number 
of EROs selected hy r — K > 5.5, (9) by i — K > 4.4, and (10) by their combination; column (11) is the surface density 
corresponding to columns (8) and (9), and (12) is the range of surface densities at if < 18 resulting from either r or j 
band based selection. 



ERO searches (eg. Small et al. I2(J(J2|I this example case 
extreme dusty ERO corresponds to a star formation rate 
of SFR ss 600 Mq year"-'^ (as calculated from Mann et al. 
120021 using total IR luminosity). At K < 18 we have 16 
IRTF EROs, and a maximum of two of them are detected 
in the mid-IR (the two ambigous cases). Since we were 
specifically targeting ISOCAM sources and thus would be 
biased to dusty EROs, we can assign a strong upper limit 
of < 13% for if < 18 EROs being strong starbursts of 
ULIRG type. 

Separation of more quiescent EROs using broad band 
colours is difficult because of the lack of optical detections 
(or J, H band data) for many of the EROs. However, we 
do note that for example EROs #15-27 in Table|21are the 
same population as that identified in Paper I as early type 
galaxies based on their red r—K but bluer i—K colours (or 
equivalently, red r — i and blue i — Z colour, see e.g. SEDs 
in Fig. 7. of Paper I). Since part of the other EROs are 
expected to also be ellipticals, the indication is totally con- 



sistent with other studies finding the early type vs. dusty 
EROs to be somewhere around 50% with a large uncer- 
tainty. We also calculated photometric redshifts for those 
of our EROs which have photometric detections in three or 
more bands using the HYPERZ code (Bolzonella, Miralles 
& Pell6 l20U0|l : These 14 EROs have redshifts ranging from 
0.7 to 1.5, with an average of z = 1.1. Nearly all of them, 
12/14, are best fitted by an evolved stellar population of 
at least 400 Myr old after an instantaneous starburst, i.e. 
essentially an elliptical. 

4.2. ERO counts and surface densities 

We detect a total ERO surface density at if < 18.5 of 
0.26—0.36 arcmin^^ in the IRTF sample, depending on the 
colour selection criteria. See Table |3| At if < 18 we have 
0.19 arcmin^^ r - if selected EROs. These IRTF counts 
are a factor of 2—3 higher than average literature counts 
using the equivalent i? — if > 5, if < 18 limits: Daddi et 
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al. H2UUUa|l give 0.08 arcmin^^^ Roche et al. 1^1717^ find 
0.05 arcmin"^ R — K > 5.0 objects. Yan & Thompson 
H2003|l give a compilation of numerous counts with differ- 
ent selection criteria, including I — K>A, I — H>3, 
and R — K > 5, which have an approximate average of 
0.1 arcmin^ at < 18. 

Some of the differences between the surface densities 
of the samples may of course be attributed to differ- 
ent selection criteria and different photometric measure- 
ment methods. Previous surveys have however found large 
field-to-field variations in ERO densities. As will be dis- 
cussed in the next section, the apparent overdensity of our 
IRTF EROs results from strong clustering of bright EROs 
around several ISOCAM targets. 



4.3. ERO excess around ISOCAM fields 

The ERO distribution is not uniform, qualitatively there 
are clearly more EROs in IRTF frames which also have 
an ISOCAM galaxy. The contrast is even larger when 
only faint MIR sources are considered. To obtain the 
result quantitavely, we calculated the surface densities 
of EROs within 90" radii around those 15//m ISOCAM 
sources which had spectroscopic or photometric redshifts 
above z > 0.5 or which did not have redshift esti- 
mates (i.e. the faintest NIR cases with no optical coun- 
terparts to WFS limits). The search radius was selected 
to roughly match sizes of central regions of rich clusters 
at z ~ 1, where 90" corresponds to ~ 0.8 Mpc. Cases 
where an ISOCAM source fell outside the IRTF frame 
but closer than 90" from our NIR data, were included 
for the appropriate part of the frame area. The nearby, 
z < 0.5 ISO-galaxies were excluded and added to the field 
sample, with one exception: we chose to assign source 
ELAISC15_J163541. 7+405900, the second "ambiguous" 
source in Sect. 14.31 for reasons discussed therein, to an 
ERO 6" away with a photometric redshift oi z — 0.85 
instead of the official ELAIS catalogue identification of a 
z = 0.188 spiral. The main results of this paper do not 
change with this decision, and results will also be given 
for the case with a contrary decision. The relevant ISO 
sources are listed in Table 01 with the corresponding num- 
ber of detected EROs. 

For the 20 ISOCAM source related EROs in an 
area of 17.2 arcmin^, we get a surface density of 
1.16 arcmin"^, whereas for the "field ERO" sample, N=20 
in 54.6 arcmin^, the surface density is 0.37 arcmin"^. 
We calculated the densities in true non-MIR frames and 
nearby-MIR frames also separately: the results are 0.35 
(N=14) and 0.38 (N=6) EROs per arcmin^, respectively. 
Changing the category of the 4 EROs found around 
ELAISC15.J163541. 7+405900, the ISO-source ERO den- 
sity would decrease to 1.06 arcmin^ and the field density 
increase to 0.43 arcmin"^. 



4.3.1. The EROs in ISO fields 

FigureQlshows the cumulative IRTF counts separated into 
the ISOCAM source related EROs (triangles) and field 
EROs (circles). 

We note that the surface density measured from the 
ISOCAM frames is fully consistent with ERO surveys per- 
formed in fields of high-z AGN, radio galaxies, or radio- 
loud QSO (eg. Wold et al. 12113 Cimatti et aL lWRTI Hall 
& GreenQ998). With this in mind, it is interesting to test 
how our EROs correlate with the ELAIS final catalogue 
radio-sources. We searched for EROs around (45" radius 
used) VLA 1.4 GHz detections with spectroscopic or pho- 
tometric redshifts of z > 0.5. The ERO density becomes 
0.80 arcmin^^, twice as high as the "field" value. We note 
that 2 of the 8 ISOCAM sources listed in Table H are 
also radio sources; also, 3 other listed cases have a radio 
source in the searched area. Thus, there clearly is no one- 
to-one corresponce of radio and mid-IR detections, though 
there is overlap. Further study of the relation of the mid- 
IR sources to the VLA sources is out of the scope of the 
present paper. 

4.3.2. The field EROs 

The derived field ERO surface density is totally consistent 
with other field ERO surveys. The surface density of EROs 
provides a constraint for passive evolution and hierarchi- 
cal galaxy formation models. However, since the present 
survey does not add significantly to the covered area of 
large random field ERO surveys at the depth of our sur- 
vey, we do not discuss the implications of surface density 
in detail here. Merely for reference, we plot a model from 
Daddi et al. i2000b|l in Fig. ^ which predicts the number 
of EROs resulting from a single e- folding time r = 0.1 
Gyr starburst at z = 3 and pure luminosity evolution 
since. It is seen that our IRTF counts, along with brighter 
counts from our Paper I and literature counts in the range 
if = 17 - 20 fit this PLE model remarkably well. 

4.4. Significance of the ERO overdensity 

Our IRTF sample is non-continuous in sky coverage and 
we are unable to perform quantitative clustering analysis 
with the sample. However, since EROs are found to be 
strongly clustered in random field surveys (e.g. Daddi et 
al. 1200 0a! Roche et al. 2002 ) , it is important to estimate 
the significance of the excess surface densities found above. 
We follow here the treatment of Wold et al. pIT)^ . The 
significance of the excess can be given by 

CTexc = (A^MIR - A^ficld)/o' (1) 

where A'mir is the total number of EROs in the fields of 
ISOCAM mid-IR galaxies, and A^ficid is the number in the 
field IRTF frames. The expected field count fluctuations 
are given by a and can be calculated from 

a' = iVfioid(l + NnaidA^C) (2) 
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ISO Object 


N EROs 


Redshift 


1 


ELAISC15_J163417.9+405653 


2 


1.0 


2 


ELAISC15_J163515. 6+405608 





2.5 


3 


ELAISC15_J163531. 1+410025 


4 


1.150 


4" 


ELAISC15.J163541. 7+405900 


4 


0.9 


5 


ELAISC15.J163543. 1+410750 


4 


— 


6 


ELAib015_J 153515. 7+404759 


3 


0.9 


7 


ELAib015_J 153555.8+405909 


2 


2.610 


8 


ELAibOi5_J 153730.9+410447 


1 




9^ 


ELAISC15_J160718.9+544406 





_ 


10 


ELAISC15_J160741. 1+550152 


1 


0.3 


11 


ELAISC15_J161013.4+550648 


1 


0.4 


12 


ELAISC15_J161101. 3+543331 





0.235 


13 


ELAISC15_J161103.0+543220 





0.147 


14 


ELAISC15_J163359.3+410917 





0.137 


15 


ELAISC15_J163414.2+410317 


4 


0.040 


16 


ELAISC15_J163629.1+411441 





0.4 


17 


ELAISC15_J163640.0+405538 





0.106 



Table 4. List of ELAIS 15 /im sources associated with our 
IRTF frames. The number of EROs within 90" from the 
mid-IR source is given. Also, we list the spectroscopic (3 
decimals) or best- fitting photometric (1 decimal) redshift 
for each source from Rowan- Robinson et al. H2004|l . The 
sources from #9 onward are those which have redshift 
z < 0.5 and were excluded from the proper "ISO-related" 
frames from which the overdensities are derived. The IRTF 
frames corresponding to all the listed sources contain an 
area of 33.1 arcmin'^. "Note: ELAIS catalogue associates 
this object with a z = 0.189 spiral, we associated it with 
a higher redshift alternative, see text. ''Note: optical data 
fall in CCD gap, thus no redshift estimate; clearly a nearby 
galaxy pair from NIR data. 
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Fig. 4. Cumulative ERO counts of the IRTF survey for 
the r — K > 5.5 selection (solid symbols): the stars show 
the total IRTF ERO counts, the triangles counts around 
faint ISOCAM sources, and the circles the "field" frame 
counts, as discussed in text. Error bars are Poissonian. 
The recent ERO counts in the fields of radio-loud QSOs 
by Wold et al. I|2()()3|l are shown as open circles and are 
seen to agree very well with our faint mid-IR related ERO 
counts. Literature field ERO counts with R — K > 5 from 
Daddi et al. H2000a|) and Roche et al. H2002|l are shown, as 
well as the i - if > 4 EROs of Yan & Thompson plH^ . 
The thick solid line is an averaged compilation of other 
field ERO counts by Wold et al. fM)^ . The dotted line 
shows for reference a best fitting PLE model (Daddi et al. 
I2000b|l discussed in Paper I. 



where A,^ is the amplitude of the field ERO angular two- 
point correlation function and A^C the "integral con- 
straint" necessary to correct for the finite survey area 
(see e.g. Daddi et al. I20()0al Roche et al. 20Q2J. We use 
the correlation amplitude A^^ values from the large area 
ERO clustering survey of Daddi et al. H2000 a'l which 
matches well our magnitude range of K < 18.5, and 
is consistent at the faint end with other determinations 
(Roche et al. "MJT, "^nUT Firth et aL l^UH^ . In addition, 
Daddi et al. (.2000al find that the following approxima- 
tion, which we also adopt, fits the actual detected clus- 
tering: C « 58 X Area~°'^, where the area is expressed in 
arcmin^ . 

Our survey area for the mid-IR frames is 17.2 arcmin^ 
and a total number of 20 EROs are found in the area, 
compared to 20 field EROs in 56.4 arcmin^: using the 
equations above we calculate a significance for the excess 
fexc = 3.0. 

To obtain a measure of quantitative uncertainty in the 
overall CTcxc value just calculated, we performed a sim- 
ple Monte Carlo bootstrap analysis (see e.g. Press et al. 
UnnS)): Each of the N=20 MIR-related EROs and the N=20 



field EROs is replaced by a randomly selected ERO from 
the original sample; the new samples are typically smaller, 
since the any given object may be picked upto N=20 times 
from the original sample. The process is repeated numer- 
ous times, and the calculation using Equations 1 and 2 is 
repeated after each round. We find a standard deviation 
of 0.59 for 10^ simulated CToxc values (the average becomes 
2.82, i.e. essentially the same as the initial value of 3.0). 
We thus conclude that the overdensity of EROs - in the ul- 
timately quite small area of 17.2 arcmin^ of the ISOCAM 
detection related IRTF frames - is statistically significant. 

Significances a^xc for other magnitude cut-offs and 
separated for r — K and i — K selections, are given in 
Table 13 The raw numbers of objects are becoming small, 
and the confidences thus lower. However, we do note that 
a stronger clustering around the ISOCAM sources is sug- 
gested for r — K EROs than iov i — K selected EROs. The 
i — K EROs have higher surface densities in the ISOCAM 
field by factors of 1.5 — 2.5 compared to field frames, 
whereas the r — K selected overdensities are by factors 
of 3 — 5. If correct, and not only result of small number 
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statistics, the effect would be consistent with the picture 
where the r — K EROs are preferentially early types, and 
thus more clustered (Daddi et al. 20D2|), and the i — K 
EROs are more biased to dusty starbursts and/or those 
galaxies with continuing star formation. 



4.5. Radial distribution 

We attempt to determine the distribution of the EROs 
making up the overdensity as a function of distance from 
the corresponding ISOCAM source. Fig. shows the sur- 
face densities of EROs calculated using the appopriate 
available areas around the ISOCAM sources. The range 
— 30", which corresponds to radial distance of 250 kpc 
at z = 1 includes 12 out of the 20 EROs. It is seen that the 
first bin, corresponding to 6 EROs closer than 15" from 
an ISOCAM galaxy, produces a very high surface density. 
Mean surface densities of the ISO related and field frames 
are also indicated. The trend is similar to that found for 
EROs by Wold et al. at slightly fainter magnitudes 

around high-z quasars, and by Best et al. H2()();-{|l for red 
R — K>4 galaxies around radio-loud AGN. We simulated 
our radial distribution again with the bootstrap method, 
by replacing the original 20 EROs, with given radial dis- 
tances, by randomly drawn EROs from the same sample: 
the error bars in Fig. [S] reflect the standard deviation of 
10^ simulated surface densities in each bin. However, we 
note that excluding object #4 in Table 01 would make the 
radial distribution flatter and significant only at ~ 2a level 
in the two first bins. 



4.6. High-z clusters of galaxies? 

We have demonstrated that fields containing faint mid-IR 
galaxies have a significant overdensity of EROs compared 
to field populations of EROs. Without spectra it is impos- 
sible to definitely show that the EROs are related to the 
mid-IR galaxy in question. Since most EROs have mag- 
nitude upper limits only, it is not possible to do colour- 
magnitude sequence analysis either. The concentrated ra- 
dial distribution nevertheless suggests a physical correla- 
tion with the EROs and the central mid-IR source. An 
alternative for a physical connection would be that the 
EROs belong to a foreground structure which enhances 
the probability of detection of background objects through 
gravitational lensing - see discussion eg. in Cimatti et al. 
iPTO l and Wold et al. ipUSjl . 

If the EROs indeed trace a cluster of galaxies, whether 
at the mid-IR source redshift or in the foreground, one 
would expect a larger fraction of EROs in a iC-selected 
sample of galaxies at that region, since clusters are known 
to have a larger fraction of ellipticals than the field (see 
e.g. Dressier Il980fl . From Tabled it is evident that this 
indeed is the case for our ISO related IRTF fields. ERO 
fractions of ^ 20% are very similar to findings of Wold et 
al. (W03] in their z - 2 QSO fields. The table also lists 
the ERO fractions in our whole survey. These values fall 
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Fig. 5. Average radial distribution of EROs in the areas 
related to an ISOCAM source. The mean field ERO sur- 
face density in our survey is plotted with the dashed line, 
and the density of ISO-frame EROs averaged over the 
whole distance range to 90" is shown with the dotted line. 
The error bars are derived from a bootstrap analysis (see 
text) - there is a significant concentration of EROs within 
30"of the ISOCAM galaxy. 



r-K >5.5 r-K> 5.8 



K > 4.4 



ISO frames 



K < 17.5 


0.09 


0.09 


0.09 


K < 18.0 


0.26 


0.19 


0.10 


K < 18.5 


0.21 


0.14 


0.11 


K < 18.7 


0.20 


0.14 


0.18 






All frames 




K < 17.5 


0.03 


0.02 


0.02 


K < 18.0 


0.09 


0.07 


0.06 


K < 18.5 


0.12 


0.07 


0.08 


K < 18.7 


0.11 


0.05 


0.10 



Table 5. ERO fractions among iiT-selected galaxies. 



consistently in between the high values of the mid-IR fields 
and those in random fields as in Daddi et al. |2000a) and 
our Paper I values, which have fractions increasing from 
0.01 to 0.08 in the tabulated magnitude range. 

5. Summary 

We have searched for EROs in a targeted near-IR survey 
in ELAIS fields, using r - K > 5.5, r - K > 5.8, and 
i — K > 4.4 colour criteria. These are equivalent to the 
commonly used R-K > 5, R-K > 5.3, and I-K > 4. We 
find 40 EROs, 14 of which are common to both r — K > 5.5 
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and i — K > 4.4 selection. We study the number counts 
and properties of these EROs, both close to and away from 
the mid-IR ELAIS sources, and summarize the results as 
follows: 

— Taking advantage of overlapping mid-IR data, we 
search for dusty EROs, since only these should be de- 
tected with the used 15fim ISOCAM band. One clear 
mid-IR ERO is found, and two other which are some- 
what ambiguous. Taking into account detection limits 
we can limit the number of strong ULIRG type star- 
bursts in the K < 18 ERO population to < 13%. 

— Broad-band colours and SED fits of those EROs with 
detections in at least 3 bands suggest these are mainly 
early type evolved galaxies at z ~ 1. 

— Cumulative number counts are provided for the EROs. 
The total counts are found to be higher than field 
counts in the literature, but this is due to excess EROs 
discovered in ISOCAM fields. When fields with no mid- 
IR sources are considered, the ERO counts fit well 
other data. Our field ERO counts, as well as other 
literature data, are very well fit by pure luminosity 
evolution models in our magnitude range. 

— We derive a significant overdensity of EROs in areas 
within 90" of faint z > 0.5 mid-IR sources. The re- 
sult is equivalent to what other studies have found in 
high redshift AGN and QSO fields. The high fraction 
of EROs among all X-selected galaxies in these fields 
suggests that the EROs are part of high redshift clus- 
ters of galaxies. The fields are small, however, and clus- 
ters and e.g. filaments of large scale structure cannot 
be properly differentiated. Furthermore, r — K selected 
EROs dX K < 18.5 are (marginally) more clustered in 
the mid-IR fields than the i — K selected EROs, con- 
sistent with a picture that the r — K selected EROs 
contain preferentially more early type galaxies than 
fhei-K EROs. 

— It is not yet clear whether the ERO overdensities are 
physically connected to the mid-IR sources (or AGN) 
or rather part of a foreground cluster or other surface 
density peak grationally enhancing the detections of 
background mid-IR sources. However, the radial distri- 
bution of EROs around the mid-IR detection suggests 
a physical connection. 
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RA (J2000) 


DEC (J2000) 




r' 


i' 


Z 


K 


1 


242.275418 


54.967102 


> 24. 


10 


> 23.20 


> 22.00 


18.38 ±0.11 


2 


242.365120 


54.469628 


> 24. 


,10 


> 23.20 


> 22.00 


18.43 ± 0.10 


3 


248.588973 


41.052093 


> 24. 


,10 


> 23.20 


21.09 ±0.09 


18.06 ± 0.08 


4 


248.593709 


41.053258 


> 24. 


,10 


> 23.20 


> 22.00 


18.07 ± 0.08 


5 


248.597116 


41.044969 


> 24. 


,10 


> 23.20 


> 22.00 


17.99 ±0.09 


6 


248.925210 


40.985268 


> 24. 


,10 


22.61 ±0.15 


21.38 ±0.12 


17.57 ±0.07 


7 


248.927333 


40.982732 


> 24. 


,10 


22.74 ± 0.16 


> 22.00 


18.25 ± 0.13 


8 


248.931334 


41.136512 


23.51 ± 0. 


,12 


22.55 ± 0.14 


> 22.00 


17.94 ± 0.08 


9 


248.964062 


41.101228 


> 24. 


,10 


> 23.20 


> 22.00 


18.02 ± 0.10 


10 


248.994868 


41.063800 


> 24. 


,10 


> 23.20 


21.57 ±0.14 


18.38 ±0.11 


11 


249.022827 


40.783118 


> 24. 


,10 


> 23.20 


> 22.00 


17.29 ± 0.06 


12 


249.059398 


40.794446 


> 24. 


,10 


> 23.20 


> 22.00 


17.92 ± 0.17 


13 


249.222717 


40.992068 


> 24. 


,10 


> 23.20 


> 22.00 


18.41 ± 0.24 


14 


249.311503 


41.158857 


> 24. 


,10 


> 23.20 


21.52 ±0.14 


18.39 ±0.12 



15 


241, 


.931449 


55 


.037148 




> 24. 


10 


22, 


,54 


± 


0. 


,10 




> 22, 


,00 


18, 


,40 


± 


0. 


,13 


16 


242, 


.382161 


54 


.471598 




> 24. 


,10 


22, 


,38 


± 


0, 


,11 


21. 


,27 ±0, 


,09 


18, 


,30 


± 


0, 


,10 


17 


242, 


.984000 


54, 


.705978 




> 24. 


,10 


22, 


.62 


± 


0, 


,13 


21, 


,13 ±0, 


.09 


18, 


.56 


± 


0, 


,15 


18 


248, 


.598040 


40, 


.953026 




> 24. 


,10 


22, 


.16 


± 


0. 


,10 


20. 


,92 ± 0, 


.08 


18, 


.08 


± 


0, 


,11 


19 


248, 


,655921 


40, 


.951751 


23, 


,20 ± 0. 


,09 


21, 


,50 


± 


0, 


,05 


20, 


,53 ±0, 


,05 


17, 


,66 


± 


0, 


,09 


20 


248, 


.870571 


41, 


.003286 




> 24. 


,10 


22, 


,07 


± 


0, 


,09 


21, 


,30 ±0, 


,10 


18, 


,52 


± 


0, 


,15 


21 


248, 


.880417 


41 


.007793 




> 24. 


,10 


21, 


,98 


± 


0, 


,09 


20, 


,69 ± 0, 


,07 


17, 


,94 


± 


0, 


,09 


22 


248, 


.880647 


41, 


.005380 




> 24. 


,10 


22, 


.32 


± 


0, 


,11 


21, 


,15 ±0, 


.09 


18, 


.11 


± 


0, 


,10 


23 


248, 


.915574 


40 


.988594 


23, 


.37 ±0. 


,11 


21, 


.63 


± 


0, 


,06 


20, 


,89 ± 0, 


.07 


17, 


.84 


± 


0, 


,08 


24 


248, 


.936320 


41, 


.151117 


23, 


,70 ± 0. 


,14 


22, 


,29 


± 


0, 


,11 




> 22, 


,00 


18, 


,17 


± 


0, 


,11 


25 


249, 


.062402 


40, 


.825955 


23, 


,19 ±0. 


,10 


21, 


,69 


± 


0, 


,07 


20, 


,77 ± 0, 


,07 


17, 


,57 


± 


0, 


,06 


26 


249, 


.311868 


41 


.159109 




> 24. 


,10 


22, 


,70 


± 


0, 


,16 




> 22, 


,00 


18, 


,36 


± 


0, 


,12 


27 


249, 


.349512 


41, 


.090600 




> 24. 


,10 


22, 


,31 


± 


0, 


,12 




> 22, 


,00 


18, 


,00 


± 


0, 


,07 



28 


242.362595 


55.183170 


23.53 ± 0. 


,14 


22.58 ± 0.09 


20.68 ± 0.07 


18.04 ± 0. 


,07 


29 


242.382952 


54.472045 


> 24. 


,10 


> 23.20 


> 22.00 


18.79 ±0. 


,13 


30 


242.570354 


55.116371 


> 24. 


,10 


23.06 ±0.14 


21.56 ±0.14 


18.64 ±0, 


,16 


31 


248.587252 


40.971250 


> 24. 


,10 


> 23.20 


> 22.00 


18.75 ±0, 


,16 


32 


248.589808 


41.050045 


> 24. 


,10 


> 23.20 


> 22.00 


18.72 ±0, 


,12 


33 


248.884767 


41.009763 


> 24. 


,10 


> 23.20 


> 22.00 


18.64 ± 0, 


,14 


34 


248.923854 


40.986371 


> 24. 


,10 


> 23.20 


> 22.00 


18.65 ± 0, 


,14 


35 


248.929511 


41.137141 


> 24. 


,10 


> 23.20 


> 22.00 


18.72 ±0, 


,13 


36 


248.930122 


41.131120 


> 24. 


,10 


> 23.20 


> 22.00 


18.73 ±0, 


,16 


37 


249.024081 


40.781571 


> 24. 


,10 


> 23.20 


> 22.00 


18.61 ±0, 


,13 


38 


249.047083 


40.814442 


> 24. 


,10 


> 23.20 


> 22.00 


18.78 ± 0. 


,14 


39 


249.103560 


41.404861 


> 24. 


,10 


> 23.20 


> 22.00 


18.79 ± 0. 


,16 


40 


249.234015 


40.996663 


> 24. 


,10 


> 23.20 


> 22.00 


18.75 ± 0. 


,32 



Table 2. Photometry of the found EROs. Sources 1-14 are those with r — K > 5.5 and i — K > 4.4, while 15-27 and 
28-40 are selected by only r — K > 5.5 and i — K < 4.4, respectively. 
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r - is: > 5.5 




r- A- > 5.8 




I ~ K > 4.4 








All frames (area 74 arcmin"^) 








K iimit 
















(mag) 


N 


arcmin"^ 


N 


arcmin"^ 


N 


• -2 

arcmm 




K < 17.5 


3 


0.041 


2 


0.027 


2 


0.027 




K < 18.0 


14 


0.189 


11 


0.149 


9 


0.122 




K < 18.5 


27 


0.365 


12 


0.162 


19 


0.257 




K < 18.7 


27 


0.365 


12 


0.162 


27 


0.365 








ISO frames (area 17.2 arcmin'^) 
















Ek 




Ek 






IN 


arcmin"^ cxexc 


N 


arcmin"^ (Texc 


IN 


arcmm 




K < 17.5 


2 


0.116 — 


1 


0.058 — 


1 


0.058 





K < 18.0 


8 


0.465 3.4 


6 


0.349 3.1 


3 


0.174 


0.6 


K < 18.5 


13 


0.756 2.6 


7 


0.407 3.8 


7 


0.407 


1.1 


< 18.7 


13 


0.756 2.9 


7 


0.407 3.7 


12 


0.698 


2.3 






field frame; 


s (area 56.4 arcmin^) 
















Eif 




Eif 






N 


arcmin"^ 


N 


arcmin"^ 


N 


. -2 

arcmm 




K < 17.5 


1 


0.018 


1 


0.018 


1 


0.018 




/s: < 18.0 


6 


0.106 


5 


0.089 


6 


0.106 




K < 18.5 


14 


0.248 


5 


0.089 


12 


0.213 




K < 18.7 


14 


0.248 


5 


0.089 


15 


0.266 





Table 3. The sample of IRTF EROs. Cumulative counts are given for different ERO colour selection criteria. 
Additionally, counts are given separately for the total survey, those areas with a faint mid-IR source within 90", 
and those without. The latter counts are intended to be taken as field ERO counts. For the ISO related counts we also 
tabulate the significance of the excess of EROs in the corresponding magnitude and colour bins. 



